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identify  and  clarify  points  of  common  misunderstanding  concerning  sampled  windows. 
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SUMMARY 


OBJECTIVE 

To  make  available  a concise  review  of  data  windows  and  their  affect  on  the  detection 
of  harmonic  signals  in  the  presence  of  broadband  noise  and  in  the  presence  of  nearby  strong 
harmonic  interference.  Also  to  call  attention  to  a number  of  common  errors  in  the  applica- 
tion of  windows  when  used  with  the  Fast  Fourier  Transform. 

RESULTS 

A comparative  list  of  common  window  performance  measures  has  been  generated.  A 
collection  of  figures  representing  classic  and  optimal  windows  is  presented.  A two-tone 
detection  experiment  is  described,  and  the  results  are  presented  in  a series  of  figures  which 
further  demonstrate  the  performance  at  the  various  windows. 

RECOMMENDATIONS 

One  of  the  optimal  or  near-optimal  windows  should  be  applied  to  data  sequences  as 
part  of  classic  harmonic  analysis.  The  optimal  windows  are,  in  fact,  families  parameterized 
over  some  index  which  allows  a degree  of  freedom  for  tailoring  a window  to  the  require- 
ments of  a given  spectral  analysis.  The  material  contained  herein  is  applicable  to  any 
processing  scheme  which  detects  and  classifies  periodic  signals  in  finite-extent  data.  In 
particular,  beamforming  is  another  area  of  high  applicability. 
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I.  INTRODUCTION 


There  is  much  signal  processing  devoted  to  detection  and  estimation.  Detection  is 
the  task  ot  determining  it  a specific  signal  set  is  present  in  an  observation,  while  estimation 
is  the  task  of  obtaining  the  values  of  the  parameters  describing  the  signal.  Often  the  signal  is 
complicated  or  is  corrupted  by  interfering  signals  or  noise.  To  facilitate  the  detection  and 
estimation  of  signal  sets,  the  observation  is  decomposed  into  a basis  set  which  spans  the  sig- 
nal space.  For  many  problems  of  engineering  interest,  the  class  of  signals  being  sought  are 
periodic  which  leads  quite  naturally  to  a spectrum  based  upon  the  simple  periodic  functions, 
the  sines  and  cosines.  Thus  the  great  theoretical  and  practical  interest  in  the  classic  Fourier 
transform. 

By  necessity,  every  observed  signal  we  process  must  be  of  finite  extent.  The  extent 
may  be  adjustable  and  selectable,  but  it  must  be  finite.  Processing  a finite-duration  observa- 
tion imposes  interesting  and  interacting  limitations  on  harmonic  analysis.  These  include, 
detectability  of  tones  in  the  presence  of  broadband  noise,  detectability  of  weak  tones  in  the 
presence  of  nearby  strong  tones,  resolvability  of  similar-strength  nearby  tones,  resolvability 
of  shifting  tones,  and  biases  in  estimating  the  parameters  of  any  of  the  aforementioned  sig- 
nals. Similar  interactions  and  limitations  apply  to  the  analysis  of  broadband  noise  signals 
such  as  those  used  in  linear  system  identification. 

For  practicality,  the  data  we  process  are  N uniformly  spaced  samples  of  the  observed 
signal.  For  convenience.  N is  highly  composite,  and  we  will  assume  N is  even.  The  harmonic 
estimates  we  obtain  through  the  discrete  Fourier  transform  (DFT)  are  N uniformly  spaced 
samples  of  the  associated  periodic  spectra.  This  approach  is  elegant  and  attractive  when  the 
processing  scheme  is  cast  as  a spectral  decomposition  in  an  N-dimensional  orthogonal  vector 
space.  The  problem  is  that  this  elegance  must  often  be  massaged  to  obtain  meaningful 
results.  We  accomplish  this  massaging  by  the  application  of  windows  to  the  sampled  data 
set  or,  equivalently,  by  smoothing  the  spectral  samples. 

The  two  operations  to  which  we  subject  the  data  are  sampling  and  windowing. 

These  can  be  performed  in  either  order.  Sampling  is  well  understood,  windowing  is  less  so, 
and  sampled  windows  for  DFTs  significantly  less  so!  We  will  address  the  interacting  con- 
siderations of  window  selection  in  harmonic  analysis  and  examine  the  special  considerations 
related  to  sampled  windows  for  DFTs. 
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II.  HARMONIC  ANALYSIS  OF  FINITE-EXTENT  DATA 
ANDTHEDFT 

Harmonic  analysis  of  finite-duration  observations  entails  the  projection  of  the 
observed  signal  on  a basis  set  spanning  the  observation  interval.  Anticipating  the  next  para- 
graph. we  define  T seconds  as  a convenient  time  interval  and  NT  seconds  as  the  observation 
interval.  The  sines  and  cosines  with  periods  equal  to  an  integer  submultiple  of  NT  seconds 
form  an  orthogonal  basis  set  for  continuous  signals  extending  over  NT  seconds.  These  are 
defined  in  Eq.  ( 1 ). 


cos  kt|  k 

[NT  J 

sin  f—  ktl  0 

[NT  J 


k = 0,l, . . . ,N-1,N,N+1 


<t<NT 


We  observe  that  by  defining  a basis  set  over  an  ordered  index  k,  we  are  defining  the  spectrum 
over  a line  (called  the  frequency  axis)  from  which  we  draw  the  concepts  of  bandwidth  and 
of  frequencies  close  to  and  far  from  a given  frequency  (which  is  related  to  resolution). 

For  sampled  signals,  the  basis  set  spanning  the  interval  of  NT  seconds  is  identical 
with  the  sequences  obtained  by  uniform  samples  of  the  corresponding  continuous  spanning 
set  up  to  the  index  N/2.  See  Eq.  (2). 


cos  knT  = cos  ^ kn  k = 0,1 , . . . ,N/2 

” " “ “ ■ 

sin  jl^knT  = sin  ^kn  n = 0,1, . . . ,N-1 


We  note  here  that  uniformly  spaced  samples  of  a continuous  orthogonal  basis  do  not,  in 
general,  form  orthogonal  sequences;  the  trigonometric  functions  (sampled  over  an  integer 
number  of  periods)  are  a convenient  exception  (rather  than  the  rule).  We  also  note  that  an 
interval  of  length  NT  seconds  is  not  the  same  as  the  interval  covered  by  N samples  separated 
by  intervals  of  T seconds.  Figure  1 demonstrates  this  by  sampling  a function  which  is  even 
about  its  midpoint  and  of  duration  NT  seconds.  The  missing  end  point  is  the  beginning  of 
the  next  period  of  the  periodic  sequence  and  is,  in  fact,  indistinguishable  from  the  zero 
point.  This  lack  of  symmetry  due  to  the  missing  (but  implied)  end  point  is  a source  of  con- 
fusion in  sampled  window  design.  This  can  be  traced  to  the  early  work  related  to  con- 
vergence factors  for  the  partial  sums  of  Fourier  Series.  The  partial  sums  (or  the  finite 
Fourier  transform)  always  include  an  odd  number  of  points  and  exhibit  even  symmetry 
about  the  origin.  Hence  the  literature  and  software  libraries  abound  with  windows  designed 
with  true  even  symmetry  rather  than  the  implied  symmetry  with  its  missing  end  point! 


FHICfDllO  FACS  BLAIK-IIOT  niKD) 


No  Sample 


58^^ 


Figure  I . N samples  of  an  even  function  taken  over  an  NT  second  interval. 

We  must  remember  for  DFT  processing  of  sampled  data  that  even  symmetry  means 
that  the  projection  upon  the  sampled  sine  sequences  is  identically  zero;  it  does  not  mean  a 
matching  left  and  right  data  point  about  the  midpoint.  To  distinguish  this  symmetry  from 
conventional  evenness  we  will  refer  to  it  as  DFT-even  (ie.,  a conventional  even  sequence  with 
the  right  end  point  removed).  Another  example  of  DFT-even  symmetry  is  presented  in 
Fig.  2 as  samples  of  a periodically  extended  triangle  wave. 

If  we  evaluate  a DFT-even  sequence  via  a finite  Fourier  transform  (by  treating  the 
+N/2  point  as  a zero-value  point),  the  resultant  continuous  periodic  function  exhibits  a non- 
zero imaginary  component.  The  DFT  of  the  same  sequence  is  a set  of  samples  of  the  finite 
Fourier  transform,  yet  these  samples  exhibit  an  imaginary  component  equal  to  zero.  Why 
the  disparity?  We  must  remember  that  the  missing  end  point  under  the  DFT  symmetry 
contributes  an  imaginary  sinusoidal  component  of  period  27t/(N/2)  to  the  finite  transform 


Periodic  extension  of 


Figure  2.  Even  sequence  under  DFT,  and  periodic  extension  of  sequence  under  DFT. 
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(corresponding  to  the  odd  component  at  sequence  position  N/2).  The  sampling  positions  of 
the  DFT  are  at  the  multiples  of  2ff/N,  which,  of  course,  correspond  to  the  zeros  of  the 
imaginary  sinusoidal  component.  An  example  of  this  fortuitous  sampling  is  shown  in  Fig.  3. 
Notice  the  sequence  f(n)  is  decomposed  into  its  even  and  odd  parts,  with  the  odd  part 
supplying  the  imaginary  sine  component  in  the  finite  transform. 


DFT-even  sequence 

1 1 

f(nl 

1 -3 

1 -2  -1 

1 0 1 

Even  component 

Ev(f(n)) 

-4  -3  -2-1  0 1 


-^n 


Odd  component 
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The  selection  of  a finite  time  interval  of  NT  seconds  and  of  the  orthogonal  trigono- 
metric basis  (continuous  or  sampled)  over  this  interval  leads  to  an  interesting  peculiarity  of 
the  spectral  expansion.  From  the  continuum  of  possible  frequencies,  only  those  which 
coincide  with  the  basis  will  project  onto  a single  basis  vector;  all  other  frequencies  will 
exhibit  non-zero  projections  on  the  entire  basis  set.  This  is  often  referred  to  as  spectral 
leakage!  Notice  this  is  a manifestation  of  processing  finite-duration  records  and  is  not  re- 
lated in  any  way  to  the  periodic  sampling. 

An  intuitive  approach  to  leakage  is  the  understanding  that  signals  with  frequencies 
other  than  those  of  the  basis  set  are  not  periodic  in  the  observation  window.  The  periodic 
extension  of  a signal  in  a period  not  commensurate  with  its  natural  period  exhibits  discon- 
tinuities at  the  boundaries  of  the  observation.  The  discontinuities  are  responsible  for 
spectral  contributions  (or  leakage)  over  the  entire  basis  set.  The  forms  of  this  discontinuity 
are  demonstrated  in  Fig.  4. 

Windows  are  applied  to  data  to  reduce  the  spectral  leakage  associated  with  finite 
observation  intervals.  From  one  viewpoint,  the  window  is  applied  to  the  data  to  reduce  the 
order  of  the  discontinuity  at  the  boundary  of  the  periodic  extension.  This  is  accomplished 
by  matching  as  many  orders  of  derivative  as  possible  at  the  boundary.  The  easiest  match,  of 
course,  is  zero.  Thus  windowed  data  are  smoothly  brought  to  zero  at  the  end  points  so  that 
the  periodic  extension  is  continuous  in  many  orders  of  derivative. 

From  another  viewpoint,  the  window  is  applied  to  the  basis  set  so  that  a signal  of 
arbitrary  frequency  will  exhibit  a significant  projection  only  on  those  basis  vectors  having  a 
frequency  dose  to  the  signal  frequency.  Of  course  both  viewpoints  lead  to  idejitical  results. 
We  can  gain  insight  into  window  design  by  occasionally  switching  between  these  two 
viewpoints. 


Figure  4.  Periodic  extension  of  sinusoid  not  periodic  in  observation  interval. 
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111.  WINDOWS  AND  FIGURES  OF  MERIT 

Windows  are  used  in  harmonic  analysis  to  deduce  the  undesirable  effects  related  to 
spectral  leakage.  Windows  impact  on  many  attributes  of  a harmonic  processor;  these  include 
detectability,  resolution,  dynamic  range,  confidence,  and  ease  of  implementation.  We  would 
like  to  identify  the  major  parameters  that  will  allow  performance  comparisons  between 
different  windows.  We  can  best  identify  these  parameters  by  examining  the  effects  on 
harmonic  analysis  of  a window. 

An  essentially  bandlimited  signal  f (t)  with  Fourier  transform  F (co)  can  be  described 
by  the  uniformly  sampled  data  set  f (nT).  This  data  set  defines  the  periodically  extended 
spectrum  pT  (co)  by  its  Fourier  series  expansion  as  identified  in  Eqs.  (3). 


F(oj)  = 


fT'(w)=  ^ f(nT)e"j‘^"T 


+1T/1 

f(t)=r  pT 

J-tIT 


(co)  e^j^^^dco/lir 


Where 


|F(co)l=0;  Icol^ 


iW 


and  where 


pT  (co)  = F (co); 


For  (real  world)  machine  processing,  the  data  must  be  of  finite  extent,  and  the  summation 
of  Eq.  (3b)  can  only  be  performed  as  a finite  sum  as  indicated  in  Eqs.  (4). 


F (u^)  ~ ^ ^ ® jwnT  . 

n=^S/2 


9 


where 


(N/2)-l 

F(tj)=  ^ f (nT)  ;Neven 

n=^/2 


(N/2)-J 


F(tj|j)=  ^ f 


(nT)  e^^^knT  ; N even 


- 27T  . 

“k-NT“ 

We  recognize  Eq.  (4a)  as  the  finite  Fourier  transform,  a summation  addressed  for  the 
convenience  of  its  even  symmetry.  Equation  (4b)  is  the  finite  Fourier  transform  with  the 
right  end  point  removed,  and  Eq.  (4c)  is  the  DFT  sampling  of  Eq.  (4b).  Of  course  for  actual 
processing,  we  require  realizability,  and  the  summation  will  shift  N/2  positions.  This  will 
affect  only  the  phase  angles  of  the  transforms,  so  for  the  convenience  of  symmetry  we  will 
address  the  windows  as  being  centered  at  the  origin.  We  also  identify  this  convenience  as  a 
major  source  of  window  misapplication.  The  shift  of  N/2  points  is  often  overlooked  or  is 
improperly  handled  in  the  definition  of  the  window.  This  is  particularly  so  when  the 
windowing  is  performed  as  a spectral  convolution.  See  the  discussion  on  the  Hanning 
window  under  the  cos®(x)  windows. 

The  question  now  posed  is,  to  what  extent  is  the  finite  summation  of  Eq.  (4a)  a 
meaningful  approximation  of  the  infinite  summation  of  Eq.  (3b)?  In  fact,  we  address  the 
question  for  a more  general  case  of  an  arbitrary  window  applied  to  the  time  function  (or 
series)  as  presented  in  Eq.  (5). 


F^(uj)  = ^ w(nT)f(nT)e"J^"T 


where 


w (nT)  = 0;  In  I > -^  , (N  even) 


w (nT)  = w (-nT);  n # ^ , w (^T)  = 0 

Let  us  now  examine  the  effects  of  the  window  on  our  spectral  estimates.  Equation 
(5)  shows  that  the  transform  (oj)  is  the  transform  of  a product.  As  indicated  in  Eq.  (6), 
this  is  equivalent  to  the  convolution  of  the  two  corresponding  transforms  (see  Appendix). 


(6) 


i 

1 


Fw(o;)  = 


(x)  W (cj-x)  dx/27T 


or 

(cj)  = F (uj)  *W  (cj) 

Equation  (6)  is  the  key  to  the  effects  of  processing  finite-extent  data.  The  equation  can  be 
interpreted  in  two  equivalent  and  enlightening  ways,  which  will  be  more  easily  visualized 
with  the  aid  of  an  example.  The  example  we  choose  is  the  sampled  rectangle  window; 
w (nT)  = 1 .0.  We  know  W (co)  is  the  Dirichlet  kernel  shown  in  Eq.  (7). 


• cjT  - I N , 
-J-7-  sinl^toT] 

W(oj)  = e ^ ? 

sin[^  0)T] 


(7) 


Except  for  the  linear  phase  shift  term  (which  will  change  due  to  the  N/2  point  shift 
for  realizability),  the  transform  has  the  form  indicated  in  Fig.  5. 


W(wl 


Figure  5.  Dirichlet  kernel  for  N point  sequence. 

The  first  observation  concerning  Eq.  (6)  is  that  the  value  of  F^  (co)  at  a particular  co,  say 
(a)  = coq.  is  the  sum  of  all  of  the  spectral  contributions  at  each  co  weighted  by  the  window 
centered  at  Uq  and  measured  at  cu.  See  Fig.  6. 
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111.  A.  EQUIVALENT  NOISE  BANDWIDTH 

From  Fig.  6 we  observe  that  the  amplitude  of  the  harmonic  estimate  at  a given  fre- 
quency is  biased  by  the  accumulated  broadband  noise  included  in  the  bandwidth  of  the 
window.  In  this  sense,  the  window  behaves  as  a filter,  gathering  contributions  for  its 
estimate  over  its  bandwidth.  For  the  harmonic  detection  problem,  we  desire  to  minimize 
this  accumulated  noise  signal,  and  we  accomplish  this  with  small-bandwidth  windows.  A 
convenient  measure  of  this  bandwidth  is  the  Equivalent  Noise  Bandwidth  (ENBW)  of  the 
window.  This  is  the  width  of  a rectangle  filter  of  the  same  peak  power  gain  that  would 
accumulate  the  same  noise  power.  See  Fig.  7. 


Figure  7.  Equivalent  noise  bandwidth  of  window. 


The  accumulated  noise  power  of  the  window  is  defined  in  Eq.  (8). 

+7t/T 


■£ 


Noise  Power  = I |W  (co)l“  dto/27T 
-tt/T 


Parseval’s  Theorem  allows  Eq.  (8)  to  be  computed  by  Eq.  (9). 


Noise  Power  = w^  (nT) 
n 


(8) 


(9) 
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Tile  peak  power  gain  of  the  window  occurs  at  to  = 0,  the  zero  frequency  power  gain,  and  is 
defined  in  Eq.  ( 1 0). 


Peak  Signal  Gain  = W (0)  w(nT) 

n 

[? 


(lOa) 


Peak  Power  Gain  = (0)  = 


(10b) 


Thus  the  equivalent  noise  bandwidth  (normalized  by  the  1 /T  bandwidth)  is  defined  in 
Eq.  (11)  and  is  tabulated  for  the  windows  of  this  report  in  Table  1 . 


(11) 

* 

111.  B.  PROCESSING  GAIN 

A concept  closely  allied  to  ENBW  is  Processing  Gain  and  Processing  loss  of  a win- 
dowed transform.  We  can  think  of  the  discrete  Fourier  transform  as  a bank  of  filters 
matched  to  the  set  of  basis  tones.  From  this  perspective  we  can  examine  the  processing  gain 
(sometimes  called  the  coherent  gain)  of  the  filter,  and  we  can  examine  the  processing  loss 
due  to  the  window  reducing  the  data  to  zero  near  the  boundaries.  Let  the  input  sampled 
sequence  be  defined  by  Eq.  (12); 

f(nT)  = Ae'^j‘^k"'^  + q (nT)  (12) 

where  q (nT)  is  a white  noise  sequence  with  variance  o^.  Then  the  signal  component  of  the 
windowed  spectrum  (the  matched  filter  output)  is  presented  in  Eq.  (13). 


F(cjk)  =y  w(nT)Ae'*'j^k"T'e"i‘^k'’’^=  aY  w(nT)  (13) 

signal  Zx 


We  see  that  the  noiseless  measurement  (the  expected  value  of  the  noisy  measurement)  is 
proportional  to  the  input  amplitude  A.  The  proportionality  factor  is  the  sum  of  the  window 
terms,  which  is  in  fact  the  DC  signal  gain  of  the  window.  For  a rectangle  window  this  factor 
is  N,  the  number  of  terms  in  the  window.  For  any  other  window,  the  gain  is  reduced  due  to 
the  window  smoothly  going  to  zero  near  the  boundaries.  This  reduction  in  proportionality 
factor  is  important  as  it  represents  a known  bias  on  spectral  amplitudes.  Coherent  power 
gain,  the  square  of  coherent  gain,  is  occasionally  the  parameter  listed  in  the  literature.  Co- 
herent gain  (the  summation  of  Eq.  (13)(  normalized  by  its  maximum  value  N is  listed  in 
Table  I. 
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The  incoherent  component  of  the  windowed  transform  is  given  in  Eq.  (14a),  and  the 
incoherent  power,  (the  mean  square  value  of  this  component)  is  given  in  Eq.  (14b). 


F(cok)  =y  w(nT)q(nT)e'j^knT 

noise  4^ 


•jiF  (c0|^) 


ss  w (nT)  w (mT)  E {q  (nT)  q*(mT)} 
n m 

= 0“^  w^  (nT) 
q L- 


Notice  the  incoherent  power  gain  is  the  sum  of  the  squares  of  the  window  terms  and  the 
coherent  power  gain  is  the  square  of  the  sum  of  the  window  terms. 

Finally,  processing  gain  (PG),  defined  as  the  ratio  of  output  signal-to-noise  ratio  to 
input  signal-to-noise  ratio,  is  given  in  Eq.  ( 1 5). 


^ w (nT)  Voq  ^ w^  (nT) 

PG  = ^o^^o  = J1 . ' n 

Si/N,  ^2|„l 


w (nT) 


^w^  (nT) 


Notice  processing  gain  is  the  reciprocal  of  the  normalized  equivalent  noise  bandwidth.  Thus 
large  ENBW  suggests  a reduced  processing  gain.  This  is  reasonable,  since  an  increased  noise 
bandwidth  permits  additional  noise  to  contribute  to  a spectral  estimate. 

in.  C.  SCALLOPING  LOSS 

An  important  consideration  related  to  minimum  detectable  signal  is  called  scalloping 
loss  or  picket  fence  effect.  We  have  considered  the  windowed  DFT  as  a bank  of  matched 
filters  and  have  examined  the  processing  gain  and  the  reduction  of  this  gain  ascribable  to  the 
window  for  tones  matched  to  the  basis  vectors.  The  basis  vectors  are  tones  with  frequencies 
equal  to  multiples  of  f^/N  (with  f^  being  the  sample  frequency).  These  frequencies  are  sam- 
ple points  from  the  spectrum  and  are  normally  referred  to  as  DFT  output  points  or  as  DFT 
bins.  We  now  address  the  question,  what  is  the  loss  in  processing  gain  for  a tone  of  fre- 
quency midway  between  two  bin  frequencies  (that  is,  at  frequencies  (k  + 1 /2)  f^/N))?  Re- 
turning to  Eq.  ( 1 3),  with  Wjj  replaced  by  w y2)>  we  determine  the  processing  gain  for 
this  half-bin  frequency  shift  as  defined  in  Eq.  (16a).  We  also  define  the  scalloping  loss  as  the 
ratio  of  coherent  gain  for  a tone  located  half  a bin  from  a DFT  sample-point  to  the  coherent 
gain  for  a tone  located  at  a DFT  sample-point,  as  indicated  in  Eq.  (16b). 
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(16a) 


t 


F (cj 


signal 


where 


w 

NT 


Scalloping  Loss 


W(0) 


(16b) 


Scalloping  loss  represents  the  maximum  reduction  in  processing  gain  due  to  signal  frequency. 
This  loss  has  been  computed  for  the  windows  of  this  report  and  has  been  included  in  Table  1 . 

111.  D.  WORST-CASE  PROCESSING  LOSS 

We  now  make  an  interesting  observation.  We  define  worst-case  processing  loss  as  the 
sum  of  maximum  scalloping  loss  of  a window  and  of  processing  loss  due  to  that  window 
(both  in  dB).  This  number  is  the  reduction  of  output  signal-to-noise  ratio  as  a result  of  win- 
dowing and  of  worst-case  frequency  location.  This  of  course  is  related  to  the  minimum 
detectable  tone  in  broadband  noise.  It  is  interesting  to  note  that  the  worst-case  loss  is 
always  between  3.0  dB  and  4.3  dB.  Windows  with  worst-case  processing  loss  exceeding 
3.8  dB  are  very  poor  windows  and  should  not  be  used.  Additional  comments  on  poor  win- 
dows will  be  found  under  Minimum- Resolution  Bandwidth  (Section  III.F).  We  can  conclude 
from  the  combined  loss  figures  of  Table  1 and  from  Fig.  1 0 that  for  the  detection  of  single 
tones  in  broadband  noise,  nearly  any  window  (other  than  the  rectangle)  is  as  good  as  any 
other.  The  difference  between  the  various  windows  is  less  than  1 .0  dB  and  for  good  win- 
dows is  less  than  0.7  dB.  The  detection  of  tones  in  the  presence  of  other  tones  is,  however, 
quite  another  problem.  Here  the  window  does  have  a marked  affect,  as  will  be  demonstrated 
shortly. 

III.  E.  SPECTRAL  LEAKAGE 


Returning  to  Eq.  (6)  and  to  Fig.  6,  we  observe  the  spectral  measurement  is  affected 
not  only  by  the  broadband  noise  spectrum,  but  also  by  the  narrowband  spectrum  which 
falls  within  the  bandwidth  of  the  window.  In  fact,  a given  spectral  component  say  at 
u)  = oJq  will  contribute  output  (or  will  be  observed)  at  another  frequency,  say  at  co  = cog 
according  to  the  gain  of  the  window  centered  at  coq  and  measured  at  cja.  This  is  the  effect 
normally  referred  to  as  spectral  leakage  and  is  demonstrated  in  Fig.  8 with  the  transform  of 
a finite-duration  tone  of  frequency  cOg 

This  leakage  causes  a bias  in  the  amplitude  and  the  position  of  a harmonic  estimate. 
Even  for  the  case  of  a single  real  harmonic  line  (not  at  a DFT  sample  point),  the  leakage 
from  the  kernel  on  the  negative  frequency  axis  biases  the  kernel  on  the  positive  frequency 
line.  This  bias  is  most  severe  and  most  bothersome  for  the  detection  of  small  signals  in  the 
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Figure  8.  Spectral  leakage  effect  of  window. 

presence  of  nearby  large  signals.  To  reduce  the  effects  of  this  bias,  the  window  should 
exhibit  low-amplitude  sidelobes  far  from  the  central  main  lobe  and  the  transition  to  the  low 
sidelobes  should  be  very  rapid.  One  indicator  of  how  well  a window  suppresses  leakage  is 
the  peak  sidelobe  level  (relative  to  the  main  Ipbe):  another  is  the  asymptotic  rate  of  falloff 
of  these  sidelobes.  These  indicators  are  listed  in  Table  1. 

III.  F.  MINIMUM  RESOLUTION  BANDWIDTH 

Figure  9 suggests  another  criterion  with  which  we  should  be  concerned  in  the  win- 
dow selection  process.  Since  the  window  imposes  an  effective  bandwidth  on  the  spectral 
line,  we  would  be  interested  in  the  minimum  separation  between  two  equal-strength  lines 
such  that  for  arbitrary  spectral  locations  their  respective  main  lobes  can  be  resolved.  The 
classic  criterion  for  this  resolution  is  the  width  of  the  window  at  the  half-power  points  (the 
3.0-dB  bandwidth).  This  criterion  reflects  the  fact  that  two  equal-strength  main  lobes 
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Non  Resolvable  Peaks  Resolvable  Peaks 

Figure  9.  Spectral  resolution  of  nearby  kernels. 


17 


separated  in  frequency  by  less  than  their  3.0-dB  bandwidths  will  exhibit  a single  spectral 
peak  and  will  not  be  resolved  as  two  distinct  lines.  The  problem  with  this  criterion  is  that  it 
does  not  work  for  the  coherent  addition  we  find  in  the  DFT.  The  DFT  output  points  are 
the  coherent  addition  of  the  spectral  components  weighted  through  the  window  at  a given 
frequency. 


If  two  kernels  are  contributing  to  the  coherent  summation,  the  sum  at  the  crossover 
point  (nominally  halfway  between  them)  must  be  smaller  than  the  individual  peaks  if  the 
two  peaks  are  to  be  resolved.  Thus,  at  the  crossover  points  of  the  kernels,  the  gain  from 
each  kernel  must  be  less  than  0.5,  or  the  crossover  points  must  occur  beyond  the  6.0  dB- 
points  of  the  windows.  Table  1 lists  the  6.0-dB  bandwidths  of  the  various  windows 
examined  in  this  report.  From  the  table,  we  see  that  the  6.0-dB  bandwidth  varies  from  1 .2 
bins  to  2.6  bins,  where  a bin  is  the  fundamental  frequency  resolution  tOg/N.  The  3.0-dB 
bandwidth  does  have  utility  as  a performance  indicator  as  shown  in  the  next  paragraph. 
Remember  however,  it  is  the  6.0-dB  bandwidth  which  defines  the  resolution  of  the  windowed 
DFT. 

From  Table  I we  see  that  the  noise  bandwidth  always  exceeds  the  3.0-dB  bandwidth. 
The  difference  between  the  two,  referenced  to  the  3.0-dB  bandwidth,  appears  to  be  a sensi- 
tive indicator  of  overall  window  performance.  We  have  observed  that  for  all  the  good  win- 
dows on  the  table,  this  indicator  was  found  to  be  in  the  range  of  4 to  5.5  percent. 

Those  windows  for  which  this  ratio  is  outside  that  range  either  have  a wide  main  lobe  or  a 
high  sidelobe  structure  and,  hence,  are  characterized  by  high  processing  loss  or  by  poor  two- 
tone  detection  capabilities.  Those  windows  for  which  this  ratio  is  inside  the  4-  to  5-percent 
range  are  found  in  the  lower  left  corner  of  the  performance  comparison  chart.  Fig.  10, 
which  is  described  next. 

While  Table  1 does  list  the  common  performance  parameters  of  the  windows  ex- 
amined in  this  report,  the  mass  of  numbers  is  not  enlightening.  We  do  realize  that  the  side- 
lobe  level  (to  reduce  bias)  and  the  worst-case  processing  loss  (to  maximize  detectability)  are 
probably  the  most  important  parameters  on  the  table.  Figure  10  shows  the  relative  position 
of  the  windows  as  a function  of  these  parameters.  Windows  residing  in  the  lower  left  corner 
of  the  figure  are  the  good  performing  windows.  They  exhibit  low  sidelobe  levels  and  low 
worst-case  processing  loss.  Of  course,  the  proof  of  the  pudding  is  in  the  eating;  see  the 
conclusion  section  and  Section  V. 
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TABLK  1 . WINDOWS  AND  I ICURKS  OI  MERIT. 


WINDOW 

HIGHEST 

SIDE- 

LOBE 

LEVEL 

(dB) 

SIDE- 

LOBE 

FALL- 

OFF 

(dB/OCT) 

~~  ■ — n 

COHERENT 

GAIN 

EQUIV. 

NOISE 

BW 

(BINS! 

3.0^B 

BW 

miNS) 

SCALLOP 

LOSS 

(dB) 

WORSE 

CASE 

PROCESS 

LOSS 

(dB) 

6.0<IB 

BW 

(BINS) 

RECTANGLE 

-13 

“6 

1.00 

' 

1.00 

0.89 

3.92 

3.92 

1.21 

TRIANGLE 

-27 

-12 

0.50 

1.33 

1.28 

1.82 

3.07 

1.78 

C0S“(X1 

a » 1.0 

-23 

-12 

0.64 

1.23 

1.20 

2.10 

3.01 

1.65 

HANNING 

a » 2.0 

-32 

-18 

0.50 

1.50 

1.44 

1.42 

3.18 

2.00 

a»3.0 

-39 

-24 

0.42 

1.73 

1.66 

1.08 

3.47 

2.32 

a»4.0 

-47 

-30 

0.38 

1.94 

1.86 

0.86 

3.75 

2.59 

HAMMING 

-43 

“6 

0.54 

1.36 

1.30 

1.78 

3.10 

1.81 

RIESZ 

“21 

-12 

0.67 

1.20 

1.16 

2.22 

3.01 

1.59 

R I EM  ANN 

-26 

-12 

0.59 

1.30 

1.26 

1.89 

3.03 

1.74 

OE  LA  VALLE 
POUSSIN 

“53 

'24 

0.38 

1.92 

1.82 

0.90 

3.72 

2.55 

TUKEY 

a = 0.25 

-14 

“18 

0.88 

1.10 

1.01 

2.96 

3.39 

1.38 

a « 0.50 

-15 

-18 

0.75 

1.22 

1.15 

2.24 

3.11 

1.57 

a > 0.75 

-19 

-18 

0.63 

1.36 

1.31 

1.73 

3.07 

1.80 

BOHMAN 

-46 

-24 

0.41 

1.79 

1.71 

1.02 

3.54 

2.38 

POISSON 

a -2.0 

-19 

-6 

0.44 

1.30 

1.21 

2.09 

3.23 

1.6B 

a -3.0 

-24 

-6 

0.32 

1.65 

1.45 

1.46 

3.64 

2.06 

a - 4.0 

-31 

-6 

0.25 

2;0e 

1.75 

1.03 

4.21 

2,58 

HANNING- 

a « 0.5 

-35 

-18 

0.43 

1.61 

1.54 

1.26 

3.33 

2.14 

POISSON 

a - 1.0 

-39 

-18 

0.38 

1.73 

1.64 

1.11 

3.50 

2.30 

a -2.0 

NONE 

-18 

0.29 

2.02 

1.87 

0.87 

3.94 

2.66 

CAUCHY 

<1-3.0 

-31 

-6 

0.42 

1.48 

1.34 

1.71 

3.40 

1.90 

a -4.0 

-35 

-6 

0.33 

1.76 

1.50 

1.36 

3.83 

2.20 

d - 5.0 

-30 

-6 

0.28 

2.06 

1.68 

1.13 

4.28 

2.53 

GAUSSIAN 

a - 2.5 

-42 

-6 

0.51 

1.39 

1.33 

1.69 

3.14 

1.86 

a - 3.0 

-55 

-6 

0.43 

1.64 

1.55 

1.25 

3.40 

2.18 

a -3.5 

-69 

-6 

0.37 

1.90 

1.79 

0.94 

3.73 

2.52 

DOLPH- 

a -2.5 

-50 

0 

0.53 

1.39 

1.33 

1.70 

3.12 

1,85 

TCHEBYSHEVa-3X) 

-60 

0 

0.48 

1.51 

1.44 

1.44 

3.23 

2.01 

a -3.5 

-70 

0 

0.45 

1.62 

1.55 

1.25 

3.35 

2.17 

a -4.0 

-80 

0 

0.42 

1.73 

1.65 

1.10 

3.48 

2.31 

KAISER- 

a -2.0 

-46 

-6 

0.49 

1.50 

1.43 

1.46 

3.20 

1.99 

BESSEL 

a -2.6 

-57 

-6 

0.44 

1.65 

1.57 

1.20 

3.38 

2.20 

a -3.0 

-69 

-6 

0.40 

1.80 

1.71 

1.02 

3.56 

2.39 

a -3.5 

-82 

-6 

0.37 

1.93 

1.83 

0.89 

3.74 

2.57 

BARCILON- 

0-3.0 

-53 

-6 

0.47 

1.56 

1.49 

1.34 

3.27 

2.07 

TEMES 

0-3.5 

-58 

-6 

0.43 

1.67 

1.59 

1.18 

3.40 

2.23 

a‘4J0 

-68 

-6 

0.41 

1.77 

1.69 

1.05 

3.52 

2.36 

HIGHEST  SIDELOBE  LEVEL,  dB 


-10 


3.0  3.5  4.0 

WORST  CASE  PROCESSING  LOSS,  dB 


Figure  10.  Comparison  of  windows;  sidelobe  levels  and  worst-case  processing  loss. 
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IV.  CLASSIC  WINDOWS 

We  will  now  catalogue  some  well-known  (and  some  not  well-known)  windows.  For 
each  window  we  will  comment  on  the  justification  for  its  use  and  identify  its  significant 
parameters.  All  the  windows  will  be  presented  as  even  (about  the  origin)  sequences  with  an 
odd  number  of  points.  To  convert  the  window  to  DFT-even,  the  right  end  point  will  be  dis- 
carded and  the  sequence  will  be  shifted  so  that  the  left  end  point  coincides  with  the  origin. 
We  will  also  use  normalized  coordinates  with  sample  period  T = 1.0,  so  that  co  is  periodic  in 
2ir  and,  hence,  will  be  identified  as  d. 

IV.  A.  RECTANGLE  WINDOW 

The  rectangle  window  is  unity  over  the  observation  interval  and  can  be  thought  of  as 
a gating  function  applied  to  the  data  so  that  they  are  of  finite  extent.  The  window  for  a 
finite  Fourier  transform  is  defined  in  Eq.  ( 1 7a)  and  is  shown  in  Fig.  1 1 . The  same  window 
for  a discrete  Fourier  transform  is  defined  in  Eq.  (17b). 


w (n)  = 1.0;  n = -y , . . . ,-1,0,1, . . 

N 

■’2 

(17a) 

w (n)  = 1.0;  n = 0,1 , . . . ,N-1 

(17b) 

The  transform  of  this  window  is  seen  to  be  the  Dirichlet  Kernel,  which  exhibits  a DFT  main- 
lobe  width  of  2(27t/N)  and  a first  sidelobe  level  approximately  13  dB  down  from  the  main- 
lobe  peak.  The  sidelobe  fall  off  at  6.0  dB  per  octave,  which  of  course  is  expected  for  a func- 
tion with  a discontinuity.  The  parameters  of  the  DFT  window  are  listed  in  Table  1 . 

With  the  rectangle  now  defined,  we  can  answer  the  question  posed  earlier:  in  what 
sense  does  the  finite  sum  of  Eq.  (1 8a)  approximate  the  infinite  sum  of  Eq.  (18b)? 

-KN/2 

F(d)=  ^ f(n)e"j"®  (18a) 

n=‘^/2 

+00 

F(0)=  ^ f(n)e-j"®  (18b) 

n=-oo 

We  observe  the  finite  sum  is  the  rectangle-windowed  version  of  the  infinite  sum.  We  recog- 
nize that  the  infinite  sum  is  the  Fourier  series  expansion  of  some  periodic  function  for 
which  the  f (n)’s  are  the  Fourier  series  coefficients.  We  also  recognize  that  the  finite  sum  is 
simply  the  partial  sum  of  the  series.  From  this  viewpoint  we  can  cast  the  question  in  terms 
of  the  convergence  properties  of  the  partial  sums  of  Fourier  series.  From  this  work  we  know 
the  partial  sum  is  the  least  mean  square  error  approximation  to  the  infinite  sum. 
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We  observe  that  mean  square  convergence  is  a convenient  analytic  concept,  but  it  is 
not  attractive  for  finite  estimates  or  for  numerical  approximations.  Mean  square  estimates 
tend  to  oscillate  about  their  mean  and  do  not  exhibit  uniform  convergence.  (The  approxi- 
mation at  a point  of  continuity  may  get  worse  if  more  terms  are  added  to  the  partial  sum.) 
We  normally  observe  this  behavior  near  points  of  discontinuity  as  the  ringing  we  call  Gibbs 
phenomenon.  It  is  this  oscillatory  behavior  we  are  trying  to  suppress  by  the  use  of  other 
windows. 

IV.  B.  TRIANGLE  (FEJER,  OR  BARTLET)  WINDOW 

The  triangle  window  for  a finite  Fourier  transform  is  defined  in  Eq.  (19a)  and  is 
shown  in  Fig.  1 2.  The  same  window  for  a DFT  is  defined  in  Eq.  ( 1 9b). 

w(n)=  1.0- -j^l;n  = -y..., -1,0,1,...,^  (19a) 

N/2  1 1 

w(n)  = ^;n  = 0,l,...,^  (19b) 

N/2  2 

= w(N-n);n  = j, ,N-1 

The  transform  of  this  window  is  seen  to  be  the  squared  Dirichlet  kernel.  Its  main-lobe 
width  is  twice  that  of  the  rectangle’s  and  the  first  sidelobe  level  is  approximately  26  dB 
down  from  the  main-lobe  peak,  again,  twice  that  of  the  rectangle’s.  The  sidelobes  fall  off 
at  -1  2 dB  per  octave,  reflecting  the  discontinuity  of  the  window  residing  in  the  first  deriva- 
tive (rather  than  in  the  function  itself).  The  triangle  is  the  simplest  window  which  exhibits 
a non-negative  transform.  This  property  can  be  realized  by  convolving  any  window  (of  half 
extent)  with  itself.  The  resultant  window’s  transform  is  the  square  of  the  original  window’s 
transform! 

A window  sequence  derived  by  self-convolving  a parent  window  contains  approxi- 
mately twice  the  number  of  samples  as  the  parent  window,  hence  about  twice  the  number 
of  degrees  of  freedom.  (Convolving  two  rectangles  each  of  N/2  points  will  result  in  a triangle 
of  N+l  points  when  the  zero  end  points  are  counted  ) The  transform  of  the  window  will 
now  exhibit  twice  as  many  zeros  as  the  parent  transform  (to  account  for  the  increased 
degrees  of  freedom).  But  how  has  the  transform  applied  these  extra  degrees  of  freedom? 

The  self-convolved  window  simply  places  repeated  zeros  at  each  location  for  which  the 
parent  transform  had  a zero.  This,  of  course,  not  only  sets  the  transform  to  zero  at  those 
points,  but  also  sets  the  first  derivative  to  zero  at  those  points.  If  the  intent  of  the  extra 
degrees  of  freedom  is  to  hold  down  the  sidelobe  levels,  then  doubling  up  on  the  zeros  is  a 
wasteful  tactic.  The  additional  zeros  might  be  better  placed  between  the  existing  zeros 
(near  the  local  peaks  of  the  sidelobes)  to  hold  down  the  sidelobes  rather  than  at  locations  for 
which  the  transform  is  already  equal  to  zero.  In  fact  we  will  observe  in  subsequent  windows 
that  very  few  good  windows  exhibit  repeated  roots. 

Backing  up  for  a moment,  it  is  interesting  to  examine  the  triangle  window  in  terms 
of  partial-sum  convergence  of  Fourier  series.  Fejer  observed  that  the  partial  sums  of  Fourier 
series  were  poor  numerical  approximations.  Fourier  coefficients  were  easy  to  generate. 
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however,  and  he  questioned  if  some  massaging  of  coefficients  might  lead  to  a new  set  with 
more  desirable  convergence  properties.  The  oscillation  of  the  partial  sum  and  the  contrac- 
tion of  those  oscillations  as  the  order  of  the  partial  sum  increased  suggested  that  an  average 
of  the  partial  sums  would  be  a smoother  function.  Figure  1 3 presents  an  expansion  of  two 
partial  sums  near  a discontinuity.  Notice  the  average  of  the  two  expansions  is  smoother 
than  either. 


Eq.  (20). 


< fN(0)  = 1 [Fn_i  (0)  + Fn_2(0)+  - . + Fo(0)]  (20) 

Where  F]^  (0)  is  the  M-term  partial  sum  of  the  series.  This  is  easily  visualized  in  Table  2, 
which  lists  the  non-zero  coefficients  of  the  first  four  partial  sums  and  their  average 
summation. 

Table  2.  Fejer  convergence  factors  as  an  average  transform. 

Fo(«)  fo 


F,  (0) 

f-1 

fo 

f+1 

F2(0) 

f-2 

f-1 

fo 

f+1 

f+2 

F3  (0) 

f-3  f-2 

f-1 

fo 

f+1 

f+2  f+3 

F^(0) 

Of 

/-4 

■f 

4C3  /-2 

3 

4‘-« 

4^0 

4f+l 

2f  K 0, 

4f+2  4f+3  4f+4 

We  see  that  the  Fejer  convergence  factors  applied  to  the  Fourier  series  coefficients  is,  in 
fact,  a Triangle  window. 
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IV.  C.  COS“(X)  WINDOWS 


This  is  actually  a family  of  windows  dependent  upon  the  parameter  a,  with  a nor- 
mally being  an  integer.  Attractions  of  this  family  include  the  ease  with  which  the  terms  can 
be  generated  and  the  easily  identified  properties  of  the  transform  of  the  cosine  function. 
These  properties  are  particularly  attractive  under  the  DFT.  The  window  for  a finite  Fourier 
transform  is  defined  in  Eq.  (21a)  and  fora  DFT  in  Eq.  (21b).  Notice  the  effect  due  to  the 
change  of  the  origin. 


w (n)  = cos  f-  ttI  ; n = - — , . . 

Ln  J 2 

..,-1,0,1,.. 

N 

■■’2 

(21a) 

w (n)  = sin  i “ 0.1.2, . 

N-1 

(21b) 

The  most  common  values  of  a are  the  integers  I through  4,  with  2 being  the  most  well 
known  (as  the  Hanning  window).  This  window  is  identified  for  values  of  a equal  to  1 and  2 
in  Eqs.  (22a),  (22b),  (23a)  and  (23b),  (the  “a”  for  the  finite  transform,  the  “b”  for  the 
DFT).  The  windows  are  shown  for  a integer  values  of  1 through  4 in  Figs.  14  through  1 7. 

a = 1 .0  (Cosine  Lobe) 

w (n)  = cos  ; n = -^, . . . 

.,-1,0,1,... 

N 

■ '2 

(22a) 

a = 1 .0  (Sine  Lobe) 

w ( n)  = sin  ; n = 0, 1 ,2, . 

. . . ,N-1 

(22b) 

a = 2.0  (Cosine  Squared,  Raised  Cosine,  Hanning) 

w (n)  = cos^ 

(23a) 

= 0.5  [i  0 + cos[^7t]] 

1 

II 

c 

,-1,0,1,.  . 

N 

' ■ '2 

a = 2.0  (Sine  Squared,  Raised  Cosine,  Hanning) 

w (n)  = sin^  ttJ 

(23b) 

'0.5[l.0-co,[|!,]] 

;n  = 0.1,2,. 

. . ,N-1 

Notice  as  a becomes  larger,  the  windows  become  smoother  and  the  transform  reflects  this 
increased  smoothness  in  decreased  sidelobe  level  and  faster  falloff  of  the  sidelobes.  but  with 
an  increased  width  of  the  main  lobe. 
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Of  parliailar  interest  in  this  family,  is  the  Hann  window  (after  the  Austrian  meteor- 
ologist, Julius  Von  Hann).*  Not  only  is  this  window  continuous,  but  so  is  its  first  derivative. 
Since  the  discontinuity  of  this  window  resides  in  the  second  derivative,  the  transform  falls 
off  as  l/cj-^  or  as  18  dB  per  octave.  Let  us  closely  examine  the  transform  of  this  window. 

We  will  gain  some  interesting  insight  and  learn  of  a clever  application  of  the  window  under 
the  DFT. 

The  sampled  Hanning  window  can  be  written  as  the  sum  of  the  sequence  as,  indicated 
in  Lq.  (24a).  Each  sequence  has  the  easily  recognized  transform  indicated  in  Eq.  (24b). 


w (n)  = 0.5  + 0.5  cos  — 
N 


.,-1,0,1,...,^ -I 


W(0)  = O.5  D(0)  + O.25 


) + D(0-h: 


where 


.6  . 

-JT  Sin 


Sin  ^0 
sin 


We  recognize  the  Dirichlet  kernel  at  the  origin  as  the  transform  of  the  constant  0.5  samples 
and  the  pair  of  translated  kernels  as  the  transform  of  the  single  cycle  of  cosine  samples. 

Note  that  the  translated  kernels  are  located  on  the  first  zeros  of  the  center  kernel  and  are 
half  the  size  of  the  center  kernel.  Also  the  sidelobes  of  the  translated  kernel  are  about  half 
the  size  and  are  of  opposite  phase  of  the  sidelobes  of  the  central  kernel.  The  summation  of 
the  three  kernels’  sidelobes,  being  in  phase  opposition,  tends  to  cancel  the  sidelobe  struc- 
ture. This  cancelling  summation  is  demonstrated  in  Fig.  1 8. 

The  partial  cancelling  of  the  sidelobe  structure  suggests  a constructive  technique  to 
define  new  windows.  The  most  well-known  of  these  is  the  Hamming  window,  which  is 
presented  in  the  next  section. 

For  the  special  case  of  the  DFT,  the  Hanning  window  is  sampled  at  multiples  of 
2rr/N,  which  of  course  are  the  locations  of  the  zeros  of  the  central  Dirichlet  kernel.  Thus 
only  three  non-zero  samples  are  taken  in  the  sampling  process.  The  positions  of  these  sam- 
ples are  at  -27r/N,  0,  and  +2jt/N.  The  value  of  the  samples  obtained  from  Eq.  (24b) 
(including  the  phase  factor  to  account  for  the  N/2  shift)  are  -1/4,  +1/2.  -1/4, 

respectively.  Note  the  minus  signs.  These  result  from  the  shift  in  the  origin  for  the  window. 
Without  the  shift,  the  phase  term  is  missing  and  the  coefficients  are  all  positive:  1/4,  1/2. 
1/4.  These  are  incorrect  for  DFT  processing,  but  they  find  their  way  into  much  of  the 
literature  and  practice. 

Rather  than  apply  the  window  as  a product  in  the  time  domain,  we  always  have  the 
option  to  apply  it  as  a convolution  in  the  frequency  domain.  The  attraction  of  the  Hanning 
window  for  this  application  is  twofold;  first,  the  window  spectra  is  non-zero  at  only  three 


*The  correct  name  of  this  window  is  “Hann.  " The  term  "Hanning"  is  used  in  this  report  to  reflect  conventional  usage. 
The  derived  term  “Hann  "d"  is  also  widely  used. 
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Figure  18.  Transform  of  Hanning  window  as  a sum  of 
three  Dirichlet  kernels. 


data  points,  and  second,  the  sample  values  are  binary  fractions,  which  can  be  implemented 
as  right  shifts.  Thus  the  Hanning-windowed  spectral  points  obtained  from  the  rectangle- 
windowed  spectral  points  are  obtained  as  indicated  in  Eq.  (25)  as  two  real  adds  and  two 
oinary  shifts  (to  multiply  by  1/2). 

Hanning  “ T ~ 2 jy  ^ ^ 

Thus  u Hanning  window  applied  to  a real  transform  of  length  N can  be  performed  as  N real 
multiplies  on  the  time  sequence  or  as  2N  real  adds  and  2N  binary  shifts  on  the  spectral  data. 
One  other  mildly  important  consideration,  if  the  window  is  to  be  applied  to  the  time  data,  is 
that  the  samples  of  the  window  must  be  stored  somewhere,  which  normally  means  addi- 
tional memory  or  hardware.  It  so  happens  that  the  samples  of  the  cosine  for  the  Hanning 
window  are  already  stored  in  the  machine  as  the  trig-table  for  the  FFT;  thus  the  window 
requires  no  additional  storage. 

IV.  D.  HAMMING  WINDOW 

The  Hamming  window  can  be  thought  of  as  a modified  Hanning  window.  (Note  the 
potential  source  of  confusion  in  the  similarities  of  the  two  names.)  Referring  back  to  Figs. 

1 5 and  1 8.  we  note  the  inexact  cancellation  of  the  sidelobes  from  the  summation  of  the 
three  kernels.  We  can  construct  a window  by  adjusting  the  relative  size  of  the  kernels  as 
indicated  in  Eq.  (26a)  to  achieve  a more  desirable  form  of  cancellation.  Perfect  cancellation 
of  the  first  .sidelobe  (at  0 = 2.5  ( 27r/N) ) corresponds  to  the  Hamming  window  as  indicated 
in  Eq.  (26b). 
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w (n)  = a + ( I-CK)  cos 


W(0)  = O(D((9)  + O.5  (1-a) 


N 


N 


w (n)  = 0.54  + 0.46  cosj—  n = — i . . . ,-1.0,1 — 

N ' 


(26a) 


(26b) 


or 


w (n)  = 0.54  - 0.46  cos 


0.1,2, ...  . N-1 


The  Hamming  window  is  shown  in  Fig.  19.  Notice  the  deep  attenuation  at  the  missing  side- 
lobe  position.  Note  also  that  the  small  discontinuity  at  the  boundary  of  the  window  has 
resulted  in  a 1/cj  (6.0  dB  per  octave)  rate  of  falloff.  The  better  sidelobe  cancellation  does 
result  in  a much  lower  initial  sidelobe  level  of -42  dB.  Table  1 lists  the  parameters  of  this 
window.  Also  note  the  loss  of  binary  weighting;  hence  the  loss  of  ease  of  a spectral  convolu- 
tion implementation. 


IV.  E.  CONSTRUCTED  WINDOWS 


Numerous  investigators  have  constructed  windows  as  products,  as  sums,  as  sections, 
or  as  convolutions  of  simple  functions  and  of  other  simple  windows.  These  windows  have 
been  constructed  for  certain  desirable  features,  not  the  least  of  which  is  the  attraction  of 
simple  functions  fe;  generating  the  window  terms.  In  general,  the  constructed  windows 
tend  not  to  be  good  windows,  and  occasionally  are  very  bad  windows.  We  have  already 
examined  some  simple  window  constructions.  The  Fejer  (Bartlett)  window,  for  instance,  is 
the  convolution  of  two  rectangle  windows;  the  Hamming  window  is  the  sum  of  a rectangle 
and  a Hanning  window;  and  the  cos'^  (X)  window  is  the  product  of  two  Hanning  windows. 
We  will  now  examine  other  constructed  windows  that  have  appeared  in  the  literature.  We 
will  present  them  so  they  are  available  for  comparison.  Later  we  will  examine  windows 
constructed  in  accord  with  some  criteria  of  optimality  (see  Sections  IV  G,  H,  I,  and  J).  Each 
window  is  identified  only  for  the  finite  Fourier  transform.  A simple  shift  of  N/2  points  and 
right  end-point  deletion  will  supply  the  DFT  version.  The  significant  figures  of  performance 
for  these  windows  are  also  found  in  Table  1. 

1 . Riesz  Window 


The  Riesz  window,  identified  in  Eq.  (27),  is  the  simplest  continuous  polynomial 
window.  It  exhibits  a discontinuous  first  derivative  at  the  boundaries;  hence  its  transform 
fall  off  like  l/cj^.  The  window  is  shown  in  Fig.  20.  The  first  sidelobe  is  -22  dB  from  the 
main  lobe.  This  window  is  similar  to  the  cosine  lobe  [Eq.  (22)]  as  can  be  demonstrated  by 
examining  its  Taylor  series  expansion. 


w (n)  = 1.0  - 


n 

N/2 


2 


0<|nl<-| 


(27) 
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2.  Riemann  Window 


I 


The  Riemann  window,  defined  in  Eq.  (28),  is  the  central  lobe  of  the  SINC  kernel. 
This  window  is  continuous,  with  a discontinuous  first  derivative  at  the  boundary.  It  is 
similar  to  the  Riesz  and  Cosine  Lobe  windows.  The  Riemann  window  is  shown  in  Fig.  2 1 . 


(28) 


3.  de  la  Valle'-Poussin  Window 

The  de  la  Valle'-Poussin  window  is  a piece-wise  cubic  curve  obtained  by  self- 
convolving  two  triangles  of  half  extent  or  four  rectangles  of  one-fourth  extent.  It  is  defined 
in  Eq.  (29). 


(29) 


The  window  is  continuous  up  to  its  third  derivative  so  that  its  sidelobes  fall  off  like  ]/u^. 
The  window  is  shown  in  Fig.  22.  Notice  the  trade-off  of  main  lobe  width  for  sidelobe  level 
Compare  this  with  the  rectangle  and  the  triangle.  It  is  a non-negative  window  by  virtue  of 
its  self-convolution  construction. 


4.  Tukey  Window 


The  Tukey  window,  often  called  the  cosine-tapered  window,  is  best  imagined  as  a 
cosine  lobe  of  width  (a/2)N  convolved  with  a rectangle  window  of  width  (1.0-  a/2)N.  Of 
course  the  resultant  transform  is  the  product  of  the  two  corresponding  transforms.  The 
window  represents  an  attempt  to  smoothly  set  the  data  to  zero  at  the  boundaries  while  not 
significantly  reducing  the  processing  gain  of  the  windowed  transform.  The  window  evolves 
from  the  rectangle  to  the  Hanning  window  as  the  parameter  a varies  from  zero  to  unity. 

The  family  of  windows  exhibits  a confusing  array  of  sidelobe  levels  arising  from  the  product 
of  the  two  component  transforms.  The  window  is  defined  in  Eq.  (30). 


w (n)  = 1.0; 


0<|n|<a- 

'J 


= 0.5 


“ 

r N 1 

n -a  T 

1 .0  + cos 

ir 

2(l-«)^ 

;«-<  |n|<- 


(30) 


The  window  is  shown  in  Figs.  23.  24,  and  25  for  values  of  a equal  to  0.25,  0.50,  and  0.75. 
respectively. 
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5.  Bohman  Window 

The  Bohmun  window  is  a construction  consisting  of  a product  of  a triangle  window 
with  a single  cycle  of  a cosine  with  the  same  period  and,  then,  a corrective  term  added  to  set 
the  first  derivative  to  zero  at  the  boundary.  Thus  the  second  derivative  is  continuous,  and 
the  discontinuity  resides  in  the  third  derivative.  The  transform  falls  off  like  \ju>^.  The  win- 
dow is  defined  in  Hq.  (3 1 ) and  is  shown  in  Fig.  26. 

“'"I"  [' 

0<|nl<-^ 

6.  Poisson  Window 

The  Poisson  window  is  a two-sided  exponential  defined  by  Eq.  (32). 

w(n)  = e N/2  0<|n(<—  (32) 

This  is  actually  a family  of  windows  parameterized  on  the  variable  a.  Since  it  exhibits  a 
discontinuity  at  the  boundaries,  the  transform  can  falloff  no  faster  than  l/w.  The  window 
is  shown  in  Figs.  27.  28.  and  29  for  values  of  a equal  to  2.0.  3.0,  and  4.0,  respectively. 
Notice  as  the  discontinuity  at  the  boundaries  becomes  smaller,  the  sidelobe  structure  merges 
into  the  asymptote.  Also  note  the  very  wide  main  lobe;  this  will  be  observed  in  Table  I as  a 
large  equivalent  noise  bandwidth  and  as  a large  worst-case  processing  loss. 

7.  Hanning-Poisson  Window 

The  Hanning-Poisson  window  is  constructed  as  the  product  of  the  Hanning  and  the 
Poisson  windows.  The  family  is  defined  in  Eq.  (33). 

0<|n|<^ 

This  window  is  similar  to  the  Poisson  window.  The  rate  of  sidelobe  falloff  is  determined  by 
the  discontinuity  in  the  first  derivative  at  the  origin  and  is  1/cu-.  Notice  as  a increases, 
forcing  more  of  the  exponential  into  the  Hanning  window,  the  zeros  of  the  sidelobe  struc- 
ture disappear  and  the  lobes  merge  into  the  asymptote.  This  window  is  shown  in  Figs.  30, 

3 1 , and  32  for  values  of  a equal  to  0.5,  1 .0,  and  2.0.  respectively.  Again  note  the  very  large 
main-lobe  width. 
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8.  Cauchy  Window 


The  Cauchy  window  is  a family  parameterized  on  a and  defined  in  Eq.  (34). 


w (n)  = 


0<|n|<y 


The  window  is  shown  in  Figs.  33,  34,  and  35  for  values  of  a equal  to  3.0,  4.0,  and  5.0, 
respectively.  Note  the  transform  of  the  Cauchy  window  is  a two-sided  exponential  (see 
Poisson  windows),  which  when  presented  on  a log  magnitude  scale  is  essentially  an  isosceles 
triangle.  This  causes  the  window  to  exhibit  a very  wide  main  lobe  and  to  have  a large 
equivalent  noise  bandwidth. 

IV.  F.  GAUSSIAN  OR  WEIERSTRASS  WINDOW 

Windows  are  smooth  positive  functions  with  tall,  thin  (ie.,  concentrated)  Fourier 
transforms.  From  the  generalized  uncertainty  principle,  we  know  we  cannot  simultaneously 
concentrate  both  a signal  and  its  Fourier  transform.  If  our  measure  of  concentration  is  the 
mean  square  time  duration  T and  the  mean  square  bandwidth  W,  we  know  all  functions 
satisfy  the  inequality  of  Eq.  (35),  with  equality  being  achieved  only  for  the  Gaussian  Pulse. 

TW^-^  (35) 

4ir 

Thus  the  Gaussian  pulse,  characterized  by  minimum  time-bandwidth  product,  is  a reasonable 
candidate  for  a window.  When  we  use  the  Gaussian  pulse  as  a window  we  have  to  truncate 
or  discard  the  tails.  By  restricting  the  pulse  to  finite  length,  the  window  no  longer  is  mini- 
mum time-bandwidth.  If  the  truncation  point  is  beyond  the  three-sigma  point,  the  error 
should  be  small  and  the  window  should  be  a good  approximation  to  minimum  time- 
bandwidth. 


The  Gaussian  window  is  defined  in  Eq.  (36). 


(n)  = e'“  [“  N/:] 


0<|nl<^ 


This  window  is  parameterized  on  a,  which  acts  as  the  reciprocal  of  the  standard  deviation,  a 
measure  of  the  width  of  its  Fourier  transform.  Increased  a will  decrease  the  width  of  the 
window  and  reduce  the  severity  of  the  discontinuity  at  the  boundaries.  This  will  result  in  an 
increased  width  transform  main  lobe  and  decreased  sidelobe  levels.  The  window  is  presented 
in  Figs.  36,  37,  and  38  for  values  of  a equal  to  2.5,  3.0,  and  3.5,  respectively.  Note  the 
rapid  drop-off  rate  of  sidelobe  level  in  the  exchange  of  sidelobe  level  for  main-lobe  width. 
The  figures  of  merit  for  this  window  are  listed  in  Table  1 . 

IV.  G.  Dolph-Tchebyshev  Window 


Following  the  reasoning  of  the  previous  section,  we  seek  a window  which,  for  a 
known  finite  duration,  in  some  sense  exhibits  a narrow  bandwidth.  We  now  take  a lead 
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from  tJie  antenna  design  people  who  have  faced  and  solved  a similar  problem.  The  problem 
is  to  illuminate  an  antenna  of  finite  aperture  to  achieve  a narrow  main-lobe  beam  pattern 
while  simultaneously  restricting  sidelobe  response.  (The  antenna  designer  calls  his  weighting 
procedure  shading.)  The  closed-form  solution  to  the  minimum  mainlobe  width  for  a given 
sidelobe  level  is  the  Dolph-Tchebyshev  window  (shading).  The  continuous  solution  to  the 
problem  exhibits  impulses  at  the  boundaries,  which  restricts  continuous  realizations  to 
approximations  (the  Taylor  approximation).  The  discrete  or  sampled  window  is  not  so 
restricted  and  the  solution  can  be  implemented  exactly. 

The  relation  T„  (X)  = COS  (nd)  describes  a mapping  between  the  nth  order 
Tchebyshev  (algebraic)  polynomial  and  the  nth  order  trigonometric  polynomial.  The  Dolph- 
Tchebyshev  window  is  defined  with  this  mapping  in  Eq.  (37)  in  terms  of  uniformly  spaced 
samples  of  the  window’s  Fourier  transform.  To  obtain  the  corresponding  window  time 
samples  w (n),  we  simply  perform  a DFT  on  the  samples  W (k)  and  then  scale  for  unity  peak 
amplitude.  The  parameter  a represents  the  log  of  the  ratio  of  main-lobe  level  to  sidelobe 
level.  Thus  a value  of  a equal  to  3.0  represents  sidelobes  3.0  decades  down  from  the  main 
lobe,  or  sidelobes  60.0  dB  below  the  main  lobe.  The  (-1  alternates  the  sign  of  successive 
transform  samples  to  reflect  the  shifted  origin  in  the  time  domain. 
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The  window  is  presented  in  Figs.  39,  40,  41,  and  42  for  values  of  a equal  to  2.5,  3.0,  3.5, 
and  4.0,  respectively.  Note  the  uniformity  of  the  sidelobe  structure;  almost  sinusoidal!  It 
is  this  uniform  oscillation  which  is  responsible  for  the  impulses  in  the  window. 

IV.  H.  KAISER-BESSEL  WINDOW 


Let  us  examine  for  a moment  the  optimality  criteria  of  the  last  two  sections.  In 
Section  F we  sought  the  function  with  minimum  time-bandwidth  product.  We  know  this  to 
be  the  Gaussian.  In  Section  G we  sought  the  function  with  restricted  time  duration,  which 
minimized  the  main-lobe  width  for  a given  sidelobe  level.  We  now  consider  a similar  prob- 
lem. For  a restricted  energy,  determine  the  function  of  restricted  time  duration  T which 
maximizes  the  energy  in  the  band  of  frequencies,  W.  Slepian,  Poliak,  and  Landau  have  deter- 
mined this  function  as  a family  parameterized  over  the  time-bandwidth  product  TW,  the 


prolate-spheroidal  wave  functions  of  order  zero.  Kaiser  has  discovered  a simple  approxima- 
tion to  these  functions  in  terms  of  the  zero-order  modified  Bessel  function  of  the  first  kind. 
The  Kaiser-Bessel  window  is  defined  in  Eq.  (38).  The  parameter  ita  is  half  of  the  time- 
bandwidth  product. 
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This  window  is  presented  in  Figs.  43,  44,  45,  and  46  for  values  of  a equal  to  2.0,  2.5,  3.0, 
and  3.5,  respectively.  Note  the  trade-off  between  sidelobe  level  and  main-lobe  width. 

IV.  1.  BARCILON-TEMES  WINDOW 

We  now  examine  the  last  criterion  of  optimality  for  a window.  We  have  already 
described  the  Slepian,  Poliak,  and  Landau  criterion.  Subject  to  the  constraints  of  fixed 
energy  and  fixed  duration,  determine  the  function  which  maximizes  the  energy  in  the  band 
of  frequencies  W.  A related  criterion,  subject  to  the  constraints  of  fixed  area  and  fixed 
duration,  is  to  determine  the  function  which  minimizes  the  energy  (or  the  weighted  energy) 
outside  the  band  of  frequencies  W.  This  is  a reasonable  criterion  since  we  recognize  that  the 
transform  of  a good  window  should  minimize  the  energy  it  gathers  from  frequencies  re- 
moved from  its  center  frequency.  Till  now,  we  have  been  responding  to  this  goal  by  maxi- 
mizing the  concentration  of  the  transform  at  its  main  lobe. 

A closed-form  solution  of  the  unweighted,  minimum-energy  criterion  has  not  been 
found.  A solution  defined  as  an  expansion  of  prolate-spheroidal  wave  functions  does  exist, 
and  it  is  of  the  form  shown  in  Eq.  (39). 


Here  the  X^n  the  eigenvalue  corresponding  to  the  associated  prolate-spheroidal  wave  func- 
tion ( \{/2n  (x,y)l  and  the  va  is  the  selected  half  time-bandwidth  product.  The  summation 
converges  quite  rapidly  and  is  often  approximated  by  the  first  term  or  by  the  first  two 
terms.  The  first  term  happens  to  be  the  solution  of  the  Slepian,  Poliak,  and  Landau  prob- 
lem, which  we  have  already  examined  as  the  Kaiser-Bessel  window. 


A closed-form  solution  of  a weighted,  minimum-energy  criterion,  presented  in  Eq. 
(40),  has  been  found  by  Barcilon  and  Temes. 
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This  criterion  is  one  which  is  a compromise  between  the  Dolph-Tchebyshev  and  the  Kaiser- 
Bessel  window  criteria. 


Like  the  DoJph-Tchebyshev  window,  the  Fourier  transform  is  more  easily  defined, 
and  the  window  time-samples  are  obtained  by  an  inverse  DFT  and  an  appropriate  scale  fac- 
tor. The  transform  samples  are  defined  in  Eq.  (41).'  See  also  Eq.  (37). 
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where 

A = sinh  (C)  = v/lO^®-  1 
B = cosh  (C) = 10® 

C = cosh-1  (10®) 

P = cosh  cj 

y(k)  = N cos'l  cos  ^ 

This  window  is  presented  in  Figs.  47,  48,  and  49  for  values  of  a equal  to  3.0,  3.5,  and  4.0, 
respectively.  The  main-lobe  structure  is  practically  indistinguishable  from  the  Kaiser- 
Bessel  main  lobe.  The  figures  of  merit  listed  on  Table  1 suggest  that  for  the  same  sidelobe 
level,  this  window  does  indeed  reside  between  the  Kaiser-Bessel  and  the  Dolph-Tchebyshev 
windows.  It  is  interesting  to  examine  Fig.  10  and  note  where  this  window  is  located  with 
respect  to  the  Kaiser-Bessel  window;  striking  similarity  in  performance! 
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Figure  23.  25%  cosine  taper  (tukey)  window, 
Fourier  transform  and  log-magnitude 
of  transform. 


Figure  24.  50%  cosine  taper  (tukey)  window, 
Fourier  transform  and  log-magnitude 
of  transform. 


Figure  35.  Cauchy  window,  Fourier  transform.  Figure  36.  Gaussian  window,  Fourier  transform, 
and  lug-magnitude  of  transform,  (a  = S.O)  log-magnitude  of  transform,  (a  = 2.5) 
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V.  HARMONIC  RESOLUTION 

We  now  describe  a simple  experiment  which  dramatically  demonstrates  the  influence 
a window  exerts  on  the  detection  of  a weak  spectral  line  in  the  presence  of  a strong  nearby 
line.  If  two  spectral  lines  reside  in  DFT  bins,  the  rectangle  window  allows  each  to  be  identi- 
fied with  no  interaction,  lb  demonstrate  this,  consider  the  signal  composed  of  two  fre-  i 

quencies  10  fg/N  and  16  tyN  (corresponding  to  the  tenth  and  the  sixteenth  DFT  bins)  and 
of  amplitudes  1.0  and  0.01  (40.0  dB  separation),  respectively.  The  power  spectrum  of  this 
signal  obtained  by  a DFT  is  shown  in  Fig.  50  as  a linear  interpolation  between  the  DFT  out- 
put points. 

We  now  modify  the  signal  slightly  so  that  the  larger  signal  resides  midway  between 
two  DFT  bins;  in  particular,  at  10.5  f^/N.  The  smaller  signal  still  resides  in  the  sixteenth  bin. 

The  power  spectrum  of  this  signal  is  shown  in  Fig.  5 1 . We  note  that  the  sidelobe  structure 
of  the  larger  signal  has  completely  swamped  the  main  lobe  of  the  smaller  signal.  In  fact,  we  l 

know  (see  Fig.  1 1 ) that  the  sidelobe  amplitude  of  the  rectangle  window  at  5.5  bins  from  ; 

the  center  is  only  25  dB  down  from  the  peak.  Thus  the  second  signal  (5.5  bins  away)  could  1 

not  be  detected  because  it  was  more  than  26  dB  down  and.  hence,  hidden  by  the  sidelobe.  j 

(The  26  dB  comes  from  the  -25-dB  sidelobe  level  minus  the  3.9-dB  processing  loss  of  the  | 

window  plus  3.0  dB  for  positive  detection.)  We  also  note  the  obvious  asymmetry  around  j 

the  main  lobe  centered  at  1 0,5  bins.  This  is  due  to  the  coherent  addition  of  the  sidelobe  | 

structures  of  the  pair  of  kernels  located  at  the  plus  and  minus  10.5  bin  positions.  We  are  1 

observing  the  self-leakage  between  the  positive  and  the  negative  frequencies.  Figure  52  is 
the  power  spectrum  of  the  signal  pair  modified  so  that  the  large-amplitude  signal  resides  at 
the  10.25-bin  position.  Note  the  change  in  asymmetry  of  the  main  lobe  and  the  reduction 
in  the  sidelobe  level.  We  still  can  not  observe  the  second  signal  located  at  bin  position  1 6.0. 

We  now  apply  different  windows  to  the  two-tone  signal  to  demonstrate  the  differ- 
ence in  second-tone  detectability.  For  some  of  the  windows,  the  poorer  resolution  occurs 
when  the  large  signal  is  at  10.0  bins  rather  than  at  10.5  bins.  We  will  always  present  the 
window  with  the  large  signal  at  the  location  corresponding  to  worst-case  resolution. 

The  first  window  we  apply  is  the  triangle  window  (see  Fig.  53).  The  sidelobes  have 
fallen  by  a factor  of  two  over  the  rectangle  windows’  lobes  (eg.,  the  -35-dB  level  has  fallen 
J to  -70  dB).  The  sidelobes  of  the  larger  signal  have  fallen  to  approximately  -43  dB  at  the 

second  signal  so  that  it  is  barely  detectable.  If  there  were  any  noise  in  the  signal,  the  second 
tone  would  probably  not  have  been  detected. 

The  next  windows  we  apply  are  the  cos®  (x)  family.  For  the  cosine  lobe,  a=  1.0, 
shown  in  Fig.  54  we  observe  a phase  cancellation  in  the  sidelobe  of  the  large  signal  located 
at  the  small  signal  position.  This  cannot  be  considered  a positive  detection.  We  also  see  the 
; spectral  leakage  of  the  main  lobe  over  the  frequency  axis.  Signals  below  this  leakage  level 

1 would  not  be  detected.  With  a = 2.0  we  have  the  Hanning  window,  which  is  presented  in 

■ Fig.  55.  We  detect  the  second  signal  and  observe  a 3.(FdB  null  between  the  two  lobes.  This 


is  still  a marginal  detection.  For  the  cos^  (x)  window  presented  in  Fig.  56,  we  detect  the 
second  signal  and  observe  a 9.0-dB  null  between  the  lobes.  We  also  see  the  improved  side- 
lobe  response.  Finally  for  the  cos^  (x)  window  presented  in  Fig.  57,  we  detect  the  second 
signal  and  observe  a 7.0-dB  null  between  the  lobes.  Here  we  witness  the  reduced  return  for 
the  trade  between  sidelobe  level  and  main-lobe  width.  In  obtaining  further  reduction  in 
sidelobe  level  we  have  caused  the  increased  main-lobe  width  to  encroach  upon  the  second 
signal. 

We  next  apply  the  Hamming  window  and  present  the  result  in  Fig.  58.  Here  we 
observe  the  second  signal  some  35  dB  down,  approximately  3.0  dB  over  the  sidelobe 
response  of  the  large  signal.  Here,  too,  we  observe  the  phase  cancellation  and  the  leakage 
between  the  positive  and  the  negative  frequency  components.  Signals  more  than  50  dB 
down  would  not  be  detected  in  the  presence  of  the  larger  signal. 

The  Riesz  window  is  the  first  of  our  constructed  windows  and  is  presented  in  Fig. 

59.  We  have  not  detected  the  second  signal,  but  we  do  observe  its  affect  as  a 20.0-dB  null 
due  to  phase  cancellation  of  a sidelobe  in  the  large  signals’  kernel. 

The  result  of  a Riemann  window  is  presented  in  Fig.  60.  Here,  too,  we  have  no 
detection  of  the  second  signal.  We  do  have  a small  null  due  to  phase  cancellation  at  the 
second  signal.  We  also  have  a large  sidelobe  response. 

The  next  window,  the  de  la  Valle'-Poussin  or  the  self-convolved  triangle,  is  shown  in 
Fig.  61.  The  second  signal  is  easily  found  and  the  power  spectrum  exhibits  a 16.0-dB  null. 
An  artifact  of  the  window  (its  lower  sidelobe)  shows  up,  however,  at  the  second  DFT  bin  as 
a signal  approximately  53.0  dB  down.  See  Fig.  22. 

The  result  of  applying  the  Tukey  family  of  windows  is  presented  in  Figs.  62,  63,  and 
64.  In  Fig.  62  (the  25-percent  taper)  we  see  the  lack  of  second-signal  detection  due  to  the 
high  sidelobe  structure  of  the  dominant  rectangle  window.  In  Fig.  63  (the  50-percent  taper) 
we  observe  a lack  of  second-signal  detection,  with  the  second  signal  actually  filling  in  one  of 
the  nulls  of  the  first  signals’  kernel.  In  Fig.  64  (the  75-percent  taper)  we  witness  a marginal 
detection  in  the  still  high  sidelobes  of  the  larger  signal.  This  is  still  an  unsatisfying  window 
because  of  the  artifacts. 

The  Bohman-construction  window  is  applied  and  presented  in  Fig.  65.  The  second 
signal  has  been  detected  and  the  null  between  the  two  lobes  is  approximately  6.0  dB.  This 
isn’t  bad,  but  we  can  still  do  better.  Note  where  the  Bohman  window  resides  in  Fig.  10. 

The  result  of  applying  the  Poisson-window  family  is  presented  in  Figs.  66,  67.  and 

68.  The  second  signal  is  not  detected  for  any  of  the  selected  parameter  values  due  to  the 
high  sidelobe  levels  of  the  larger  signal.  We  anticipated  this  poor  performance  in  Table  1 by 
the  large  difference  between  the  3.0-dB  and  the  equivalent  noise  bandwidths. 

The  result  of  applying  the  Hanning-Poisson  family  of  windows  is  presented  in  Figs. 

69,  70,  and  71.  Here,  too,  the  second  signal  is  either  not  detected  in  the  presence  of  the 
high  sidelobe  structure  or  the  detection  is  bewildered  by  the  artifacts. 
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The  Cauchy-famiJy  windows  have  been  applied  and  the  results  are  presented  in  Figs. 
72,  73,  and  74.  Here  too  we  have  a lack  of  satisfactory  detection  of  the  second  signal  and 
the  poor  sidelobe  response.  This  was  predicted  by  the  large  difference  between  the  3.0  dB 
and  the  equivalent  noise  bandwidths  as  listed  in  Table  1. 

i 

We  now  apply  the  Gaussian  tamily  ot  windows  and  present  the  results  in  Figs.  75,  76, 
and  77.  The  second  signal  is  detected  in  all  three  figures.  We  note  as  we  further  depress  the 
sidelobe  structure  to  enhance  second-signal  detection,  the  null  deepens  to  approximately 

16.0  dB  and  then  becomes  poorer  as  the  main-lobe  width  increases  and  starts  to  overlap  the 
lobe  of  the  smaller  signal. 

The  Dolph-Tchebyshev  family  of  windows  is  presented  in  Figs.  78  through  82.  We 
observe  positive  detection  of  the  second  signal  in  all  cases,  but  it  is  distressing  to  see  the 
uniformly  high  sidelobe  structure.  Here,  we  again  see  the  coherent  addition  of  the  sidelobes 
I from  the  positive  and  negative  frequency  kernels.  Notice  that  the  smaller  signal  is  not  40  dB 

down  now.  What  we  are  seeing  is  the  scalloping  loss  of  the  large  signals’  main  lobe  being 
sampled  off  of  the  peak  and  being  referenced  as  zero  dB.  Figures  78  and  79  demonstrate 
the  sensitivity  of  the  sidelobe  coherent  addition  to  main-lobe  position.  In  Fig.  78  the  larger 
signal  is  at  bin  10.5;  in  Fig.  79  it  is  at  bin  10.0.  Note  the  difference  in  phase  cancellation 
near  the  base  of  the  large  signal.  Figure  81,  the  70-dB  sidelobe  window,  exhibits  an  18-dB 
null  between  the  two  main  lobes  but  the  sidelobes  have  added  constructively  (along  with 
the  scalloping  loss)  to  the  -62.0-dB  level.  In  Fig.  82  we  see  the  80-dB  sidelobe  window 
exhibited  sidelobes  below  the  70-dB  level  and  still  managed  to  hold  the  null  between  the 
two  lobes  to  approximately  1 8.0  dB. 

j The  Kaiser-Bessel  family  is  presented  in  Figs.  83  through  86.  Here,  too,  we  have 

i positive  second-signal  detection.  Again,  we  see  the  effect  of  trading  increased  main-lobe 

width  for  decreased  sidelobe  level.  The  null  between  the  two  lobes  reaches  a maximum  of 

22.0  dB  as  the  sidelobe  structure  falls  and  then  becomes  poorer  with  further  sidelobe  level 
improvement.  Note  that  this  window  can  maintain  a 20.0-dB  null  between  the  two  signal 
lobes  and  still  hold  the  leakage  to  more  than  70  dB  down  over  the  entire  spectrum. 

Figures  87,  88,  and  89  present  the  performance  of  the  Barcilon-Temes  window. 

Note  the  positive  detection  of  the  second  signal.  There  are  slight  sidelobe  artifacts.  The 
window  can  maintain  a 20.0-dB  null  between  the  two  signal  lobes.  The  performance  of  this 
window  is  slightly  shy  of  that  of  the  Kaiser-Bessel  window,  but  the  two  are  remarkably 
similar. 
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Figure  87.  Barcilon-Temes  window  (o  = 3 .0). 
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Figure  88.  BarcUon-Temes  window  (a  = 3,5). 


Figure  89.  Barcilon-Temes  window  (a  * 4.0). 


J 


64 


VI.  CONCLUSIONS 


We  have  examined  some  classic  windows  and  some  windows  which  satisfy  some 
criteria  of  optimality.  In  particular,  we  have  described  their  effects  on  the  problem  of 
general  harmonic  analysis  of  tones  in  broadband  noise  and  of  tones  in  the  presence  of  other 
tones.  We  have  observed  that  when  the  DFT  is  used  as  a harmonic  energy-detector,  the 
worst-case  processing  loss  due  to  the  windows  appears  to  be  lower  bounded  by  3.0  dB  and 
(for  good  windows)  upper  bounded  near  3.75  dB.  This  suggests  that  the  choice  of  particular 
windows  has  very  little  effect  on  worst-case  performance  in  DFT  energy  detection.  We  have 
concluded  that  a good  performance  indicator  for  the  window  is  the  difference  between  the 
equivalent  noise  bandwidtii  and  the  3.0-dB  bandwidth  normalized  by  the  3.0-dB  bandwidth. 
The  windows  which  perform  well  (as  indicated  in  Fig.  10)  exhibit  values  of  this  ratio  be- 
tween 4.0  and  5.5  percent.  The  range  of  this  ratio  for  the  windows  listed  in  Table  1 is  3.2 
to  22.9  percent. 

For  multiple-tone  detection  via  the  DFT.  the  window  employed  does  have  a con- 
siderable effect.  Maximum  dynamic  range  of  multi-tone  detection  requires  the  transform  of 
the  window  to  exhibit  a highly  concentrated  central  lobe  with  very  low  sidelobe  structure. 
We  have  demonstrated  that  many  classic  windows  satisfy  this  criterion  with  varying  degrees 
of  Success  and  some  not  at  all.  We  have  demonstrated  the  optimal  windows  perform  best  in 
positive  detection  of  nearby  tones  of  significantly  different  amplitudes.  Also  for  the  same 
dynamic  range,  the  three  optimal  windows  (Kaiser-Bessel,  Dolph-Tchebyshev,  and  Barcilon- 
Tetnes)  are  roughly  equivalent  with  the  Kaiser-Bessel,  demonstrating  minor  performance 
adlantages  over  the  other  two.  We  note  that  while  the  Dolph-Tchebyshev  window  appears 
to  be  the  best  window  by  virtue  of  its  relative  position  in  Fig.  10,  the  coherent  addition  of 
its  constant-level  sidelobes  detracts  from  its  performance  in  multi-tone  detection.  Also  the 
sidelobe  structure  of  the  Dolph-Tchebyshev  window  exhibits  extreme  sensitivity  to  co- 
efficient errors.  This  would  affect  its  performance  in  machines  operating  with  single- 
precision fixed-point  arithmetic.  This  suggests  that  the  Kaiser-Bessel  or  the  Barcilon-Temes 
window  be  declared  the  top  performer.  My  preference  is  the  Kaiser-Bessel  window.  Among 
other  reasons,  the  coefficients  are  easy  to  generate  and  the  trade-off  of  sidelobe  level  as  a 
function  of  time-bandwidth  product  is  fairly  simple. 

We  have  called  attention  to  a persistent  error  in  the  application  of  windows  when 
performing  convolution  in  the  frequency  domain,  i.e.,  the  omission  of  alternating  signs  on 
the  window  sample  spectrum  to  account  for  the  shifted  time  origin.  We  have  also  identified 
and  clarified  a source  of  confusion  concerning  the  evenness  of  windows  under  the  DFT. 


Finally,  we  comment  that  all  of  the  conclusions  presented  about  the  window 
performance  in  spectral  analysis  are  also  applicable  to  shading  for  array  processing  of  spatial 
sampled  data,  including  FFT  beamforming. 
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APPENDIX 

THE  EQUIVALENCE  OF  WINDOWING  IN  THE  TIME  DOMAIN  TO 
CONVOLUTION  IN  THE  FREQUENCY  DOMAIN 


f(t)=  / F (co)  dw/27T 


+N/2 

W (cj)  = ^ w (nT) 

n=-N/2 


-f-WW 

F^(w)=  2 w(nT)f(nT)e+j‘^"T 


becomes 


+00  r+oo 

X / F (X)  e'jxnT  e+j‘^"T 

n=— oo 


/+«»  +00 

F(x)  2 

> n=— oo 


w (nT)  dx/27r 


-+00  +W2 

F (X)  ^ w (nT)  dx/27r 


n=-N/2 


F (x)  W (cj-x)  dx/27T 


F„  (cj)  = F(cj)  *W(w) 
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